Dynamical charge inhomogeneity and crystal-field fluctuations for 4f ions in high-T c 

cuprates 
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The main relaxation mechanism of crystal-field excitations in rare-earth ions in cuprates is be- 
lieved to be provided by the fluctuations of crystalline electric field induced by a dynamic charge 
inhomogeneity generic for the doped cuprates. We address the generalized granular model as one 
of the model scenario for such an ingomogeneity where the cuprate charge subsystem remind that 
of Wigner crystal with the melting transition and phonon-like positional excitation modes. Formal 
description of R-ion relaxation coincides with that of recently suggested magnetoelastic mechanism. 
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Inelastic neutron scattering (INS) spectroscopy is a 
powerful tool to determine unambiguously the Stark mul- 
tiplet structure and crystal-field (CF) potential in rare- 
earth (R) based high-T c superconducting materials such 
as Yi_ 2 R 2; Ba2Cu306+y. 1 ' 2 This technique provides de- 
tailed information on the electronic ground state of the 
R-ions which is important to understand the thermody- 
namic magnetic properties as well as the observed coex- 
istence between superconductivity and long-range mag- 
netic ordering of the R-ion sublattice at low tempera- 
tures. Moreover, INS spectroscopy is addressed to be 
a powerful tool for a quantitative monitoring the decay 
of the antiferromagnetic state of the parent compound 
as well as the evolution of the superconducting state 
upon doping, since the linewidths of CF transitions are 
believed to directly probe the electronic susceptibility. 
The relaxation behavior appears to be extremely depen- 
dent upon the energy at which the susceptibility is being 
probed. The crystal-field INS spectroscopy is widely used 
to reveal the opening of an electronic gap in the normal 
state of underdoped superconductors 1 - and examine its 
anisotropy 3,4 Recently, the Ho 3+ CF-INS spectroscopy 
was used to investigate the oxygen and copper isotope 
effects on the pseudogap in the high-temperature super- 
conductors Ho-124 and (LaHoSr^CuC^i^ However, the 
mechanism of the relaxation of rare-earth ions in cuprates 
becomes the issue of hot debates*^ that questions the 
current interpretation of information detected by INS 
spectroscopy. 

In the normal state the excited crystal-field levels of R- 
ion interact with phonons, spin fluctuations, and charge 
carriers. These interactions limit the life-time of the ex- 
citation; thus the observed crystal-field transitions ex- 
hibit line broadening. Similarly the conventional Fermi- 
liquid metals the interaction with the charge carriers is 
by far considered to be the dominating relaxation mech- 
anism in cuprates. This interaction is usually assumed 
to be isotropic exchange coupling with effective spin- 
Hamiltonian H ex — —2I(gj — l)(s ■ J), where I is ex- 
change integral that should be nearly independent of the 
particular R-ion under consideration, gj Lande factor, 
s is the spin momentum of a charge carrier, and J the 
total momentum of R-ion. Such a scenario seems to be 
a rather natural, if taking into account the predominant 



spin channel of neutron scattering. The detailed theory 
of the respective relaxation mechanism was developed by 
Becker, Fulde and Keller (BFK-model)i^The correspond- 
ing intrinsic linewidth appears to increase almost linearly 
with temperature (T(T) ex p 2 T) according to the well- 
known Korringa law^fi Here p is the coupling constant, 
p = I(gj - l)N(E F ), where N(E F ) is the density of 
states (DOS) at the Fermi level. Namely the deviation 
from a linear temperature dependence at low tempera- 
tures has been usually interpreted in terms of the open- 
ing of a (pseudo)gap and the associated reduction in the 
damping. Fitting the high-temperature linewidth data 
in frames of simple or modified Korringa law one obtains 
the values of coupling constant which typically vary from 
0.003 to O.OOe, 1 ! 3 ! 4 ! 5 ! 6 

It should be emphasized that the spin channel of re- 
laxation implies directly the relevance of the Fermi liquid 
scenario for cuprates ignoring many signatures of non- 
Fermi liquid behavior. However, it is unlikely to be a 
shortcoming of either model approach whether it were 
intrinsically self-consistent. However, the spin-exchange 
model has a number of visible inconsistencies, firstly 
as concerns the magnitude of coupling constant. In- 
deed, a linear temperature dependence of the relaxation 
time above T c observed in EPR studies of S-ion Gd 3+ 
in YBa2Cu307 after Korringa fitting yields the magni- 
tude of exchange integral I « 3 x 10~ 4 eV (Ref ill) 
that directly points to an unrealistically big values of 
spin coupling constants p found in all the INS experi- 
ments on CF transitions. Some problems exist with a 
Lande factor oc (gj — 1) scaling. So, Mukherjee et al~ 
when studying the system Yi_ a; R a ;Ba2Cu306+j, (R=Er, 
Ho, Tm) found | p ^ "j 1 ~ 2, instead of theoretically ex- 



I p(Tm) | 



4.5, instead of ex- 



p(Er) 

This clear disagreement evidences 
The spin-exchange sce- 



pected (g^-i) = I' and 
pected = 
against exchange mechanism, 
nario fails to explain the "strange" doping dependence 
of Tm 3+ relaxation in Tm-123i 3 « and Nd 3+ relaxation in 
(LaSrNd) 2 Cu0 4 ^ 

Finally, Staub et. aliU have found that the Lorentzian 
linewidth of the quasi-clastic neutron scattering for Tb 3+ 
in YBa2Cu307 can be properly described by a simple 
(exp(A/kBT) — 1) _1 law typical for Orbach processes 
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governed by lattice vibrations. They have shown that 
such an interpretation also describes the results obtained 
earlier on Ho 3+ and Tm 3+ . They conclude that the 
interactions with the charge carriers are negligible and 
that the interactions with the lattice vibrations are re- 
sponsible for the relaxation behaviour of the 4f electrons 
in cuprates. Therefore, the INS results which claim to 
probe the superconducting gap or the pseudo-gap should 
be re-examined in terms of Orbach processes. Similar 
conclusion was drawn by Roepke et. alm^ for Nd 3+ re- 
laxation in (LaSrNd)2Cu04. Lovesey and Staubi& have 
shown that the dynamic properties of the lanthanidc ions 
(Tb 3+ , Ho 3+ , and Tm 3+ are adequately described by a 
simple three-state model not unlike the one introduced 
by Orbach for the interpretation of electron paramag- 
netic resonance signals from a lanthanide ion in dilute 
concentration in a salt. The cross section for inelastic 
scattering of neutrons by the lanthanide ion is derived 
by constructing a pseudospin S=l model and treating 
the magnetoelastic interaction as a perturbation on the 
three crystal-field states. The scattering of neutrons is 
thus a quasielastic process and the relaxation rate is pro- 
portional to (exp(A/kBT) — 1) , where A is the energy 
of the intermediate crystal- field state at which the density 
of phonon states probed. However, this very attractive 
scenario also faces some visible difficulties with the expla- 
nation, for instance, of the unusual nonmonotonic tem- 
perature dependences and too large oxygen isotope effect 
in the INS spectra of Ho-124 and Ho-214 systems?^ some 
doping dependences in Nd-214 system^ 

On comparing two mechanisms we should underline 
their difference which seems to be of primary impor- 
tance: the spin-channel mechanism takes into account 
the fluctuations of effective magnetic field on R-ion, 
while the phonon (magnetoelastic) mechanism deals with 
that of electric field. Moreover, the conventional spin- 
channel mechanism actually probes spin fluctuations 
rather than charge fluctuations, albeit its contribution to 
the linewidth T(T) cx (I N(Ep) 2 is believed to strongly 
depend on the density of carriers. However, this rela- 
tionship is derived in frames of Fermi liquid scenario, 
and would be modified, if one addresses the typical an- 
tiferromagnetic insulating state. Interestingly, that in 
Refsl Ml^lllfil the phonon (magnetoelastic) mechanism is 
addressed as an alternative to the charge fluctuations. 
As an example, the authors point to insulating materials 
where "...the density of carriers is essentially zero...", 8 
that forbids the charge fluctuation channel of relaxation. 

We would like to emphasize the fact that both groups 
of researchers underestimate the role of the R-ion re- 
laxation due to a conventional spinless charge fluctua- 
tion channel. Indeed, the CF Hamiltonian for R-ion in 
cuprate can be written in its standard form as follows: 

k=2Afi -k<q<k 

Here 0\ are Stevens equivalent operators, Bk q = 
bkq{r k )"fk, where bk q are CF parameters, 72 = a,7 4 = [3, 



76 = 7 {a, 13, 7 are Stevens parameters); 

bkq = (bkq) + A6 fcg , 

which may be expressed in frames of the well-known 
point-charge model as follows: 

where C k is the spherical harmonics, hi(t) the charge 
number operator. Conventional metals are character- 
ized by a very short-time charge dynamics that allows 
to neglect the contribution of charge fluctuations to the 
relaxation of R-ions in the low-energy range of CF ener- 
gies, and deal with a mean homogeneous charge distri- 
bution. An altogether different picture emerges in the 
case of cuprates where we deal with various manifesta- 
tions of static and dynamic charge inhomogeneity(see, 
e.g. Refs ll^lT^ and references herein). Moreover, the INS 
spectroscopy of CF excitations itself yields the impressive 
picture of charge inhomogeneity in 123 systemji^ where 
it was found that the observed CF spectra separate into 
different local components whose spectral weights dis- 
tinctly depend on the doping level, i.e., there is clear 
experimental evidence for cluster formation. The onset 
of superconductivity can be shown to result from perco- 
lation which means that the superconductivity is an in- 
homogeneous materials property. It seems probable that 
the dynamical rearrangement of the charge system at the 
temperatures above T c somehow affects the R-ion relax- 
ation. 

At present the stripe model of inhomogeneityi^ be- 
came the most popular. This model is based on the more 
universal idea of topological phase separation, when the 
doped particles are assumed to localize inside the domain 
walls of a bare phase. 

Below we address one of the topological phase sepa- 
ration scenario which may be termed as a generalized 
granular model for doped cuprates. We assume that the 
Cu02 layers in parent cuprates may gradually loose its 
stability under electron/hole doping, while a novel self- 
organized multigranular 2D phase is believed to become 
stable. Such a situation resembles in part that of granu- 
lar superconductivity. 

The most probable possibility is that every microgram 
accumulates one or two particles. Then, the number of 
such entities in a multigranular texture nucleated with 
doping has to nearly linearly depend on the doping. Gen- 
erally speaking, each individual micrograin may be char- 
acterized by its position, nanoscale size, and the orienta- 
tion of U(l) degree of freedom. In contrast with the uni- 
form states the phase of the superfluid order parameter 
for micrograin is assumed to be unordered. The granular 
structure must be considered as being largely dynamic in 
nature. 

In the long-wavelength limit the off-diagonal ordering 
can be described by an effective Hamiltonian in terms 
of U(l) (phase) degree of freedom associated with each 
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micrograin. Such a Hamiltonian contains a repulsive, 
long-range Coulomb part and a short-range contribution 
related to the phase degree of freedom. The latter term 
can be written out in the standard for the XY model 
form of a so-called Josephson coupling 

Hj = -^2 J n cos (^ ~ <ft); (!) 

where ipi , ipj are global phases for micrograins centered 
at points respectively, Jij Josephson coupling pa- 
rameter. Namely the Josephson coupling gives rise to 
the long-range ordering of the phase of the superfluid 
order parameter in such a multi-center texture. Such a 
Hamiltonian represents a starting point for the analysis 
of disordered superconductors, granular superconductiv- 
ity, insulator-superconductor transition with {i, j) array 
of superconducting islands with phases tpi , ipj . 

To account for Coulomb interaction and allow for 
quantum corrections we should introduce into effective 
Hamiltonian the charging energy^ 

where ni is a number operator for particles bound in i- 
th micrograin; it is quantum-mechanically conjugated to 
ip: rii = —id/d(fi, (C -1 )^ stands for the capacitance 
matrix, q for a particle charge. 

Such a system appears to reveal a tremendously rich 
quantum-critical structure. In the absence of disorder, 
the T = phase diagram of the multigranular system im- 
plies either triangular, or square crystalline arrangements 
with possible melting transition to a liquid. It should be 
noted that analogy with charged 2D Coulomb gas implies 
the Wigner crystallization of multigranular system with 
Wigner crystal (WC) to Wigner liquid melting transi- 
tion, respectively. Naturally, that the additional degrees 
of freedom for micrograin provide a richer physics of such 
lattices. For a system to be an insulator, disorder is 
required, which pins the multigranular system and also 
causes the crystalline order to have a finite correlation 
length. Traditional approach to a Wigner crystallization 
implies the formation of a WC for densities lower than a 
critical density, when the Coulomb energy dominates over 
the kinetic energy. The effect of quantum fluctuations 
leads to a (quantum) melting of the solid at high densi- 
ties, or at a critical lattice spacing. The critical properties 
of a two-dimensional lattice without any internal degree 
of freedom are successfully described applying the BKT 
theory to dislocations and disclinations of the lattice, and 
proceeds in two steps. The first implies the transition 
to a liquid-crystal phase with short-range translational 
order, the second does the transition to isotropic liquid. 
For such a system provided the micrograin positions fixed 
at all temperatures, the long-wave-length physics would 
be described by an antiferromagnetic XY model with 
expectable BKT transition and gapless XY spin-wave 
mode. 



The low temperature physics in a multigranular sys- 
tem is governed by an interplay of two BKT transitions, 
for the U(l) phase and positional degrees of freedom, 
respectively.^ 2 ! Dislocations lead to a mismatch in the 
U(l) degree of freedom, which makes the dislocations 
bind fractional vortices and lead to a coupling of transla- 
tional and phase excitations. Both BKT temperatures ei- 
ther coincide (square lattice) or the melting one is higher 
(triangular lattice) i2i 

Quantum fluctuations can substantially affect these re- 
sults. Quantum melting can destroy U(l) order at suf- 
ficiently low densities where the Josephson coupling be- 
comes exponentially small. Similar situation is expected 
to take place in the vicinity of structural transitions in 
a multigranular crystal. With increasing the micrograin 
density the quantum effects result in a significant lower- 
ing of the melting temperature as compared with classical 
square-root dependence. The resulting melting temper- 
ature can reveal an oscilating behavior as a function of 
particle density with zeros at the critical (magic) densi- 
ties associated with structural phase transitions. 

In terms of our model, the positional order corresponds 
to an incommensurate charge density wave, while the 
U(l) order does to a superconductivity. In other words, 
we arrive at a subtle interplay between two orders. The 
superconducting state evolves from a charge order with 
Tc < T m , where T m is the temperature of a melting 
transition which could be termed as a temperature of 
the opening of the insulating gap (pseudo-gap!?). 

The normal modes of a dilute multigranular system in- 
clude the pseudo-spin waves propagating in-between the 
micrograins, the positional fluctuations, or quasi-phonon 
modes, which are gapless in a pure system, but gapped 
when the lattice is pinned, and, finally, fluctuations in 
the U(l) order parameter. 

The orientational fluctuations of the multigranular sys- 
tem are governed by the gapless XY models The rel- 
evant model description is most familiar as an effective 
theory of the Josephson junction array. An important 
feature of the model is that it displays a quantum-critical 
point. 

The low-energy collective excitations of multigranular 
liquid includes an usual longitudinal acoustic phonon- 
like branch. The liquid crystal phases differ from the 
isotropic liquid in that they have massive topological ex- 
citations, i.e., the disclinations. One should note that 
the liquids do not support transverse modes, these could 
survive in a liquid state only as overdamped modes. So 
that it is reasonable to assume that solidification of the 
skyrmion lattice would be accompanied by a stabiliza- 
tion of transverse phonon-like modes with its sharpening 
below melting transition. In other words an instability of 
transverse phonon-like modes signals the onset of melt- 
ing. The phonon-like modes in skyrmion crystal have 
much in common with usual phonon modes, however, 
due to electronic nature these can hardly be detected if 
any by inelastic neutron scattering. 

A generic property of the positionally ordered 



4 



skyrmion configuration is the sliding mode which is usu- 
ally pinned by the disorder. The dcpinning of sliding 
mode(s) can be detected in a low-frequency and low- 
temperature optical response. 

It should be noted that as regards the CF fluctua- 
tions, there is no principal difference between the contri- 
butions of real phonon modes and quasi-phonon modes 
of a multigranular system. Moreover, it is worth noting 
that the charge inhomogeneity in a multigranular sys- 
tem is proned to be closely coupled with lattice struc- 
tural distorsions. However, the stabilization of transverse 
phonon-like modes in multigranular system that accom- 
panies its solidification at the temperatures above T c may 
strongly affect the CF relaxation due to a mechanism 
identical to that of proposed by Lovesey and Staub (mag- 



netoelastic mechanism). In a sense, such a conclusion rec- 
onciles "old" spin-fluctuationi^ and "new" phonor*i£ii& 
approaches to the INS spectroscopy of cuprates with R- 
ions. 

Concluding we argue that the crystal-field fluctuations 
induced by the dynamic charge inhomogeneity in copper- 
oxygen subsystem may be one of the main origin of the 
broadening of the linewidth of CF transitions for 4f ions 
in high-T c cuprates providing the detection of a charge 
rearrangement accompanying the approach to T c . 
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